Generation of harmonics by a focused laser beam in the vacuum 
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We consider generation of odd harmonics by a super strong focused laser beam in the vacuum. The 
process occurs due to the plural light-by- light scattering effect. In the leading order of perturbation 
theory, generation of (2fe -|- l)th harmonic is described by a loop diagram with {2k + 2) external 
incoming, and two outgoing legs. A frequency of the beam is assumed to be much smaller than the 
Compton frequency, so that the approximation of a constant uniform electromagnetic field is valid 
locally. Analytical expressions for angular distribution of generated photons, as well as for their 
total emission rate are obtained in the leading order of perturbation theory. Influence of higher- 
order diagrams is studied numerically using the formalism of Intense Field QED. It is shown that 
the process may become observable for the beam intensity of the order of lO^^W/cm . 

PACS numbers: 12.20.Ds,11.80.La,41.85.-p,42.65.Ky 



In a recent paper [l| T. Tajima and G. Mourou sug- 
gested a path to reach an extremely high intensity level 
10^^~^*W/cm^ already in the coming decade, taking ad- 
vantage of the megajoule facilities, see also Refs. |3,b|- 
The field strength for such lasers will be very close to 
the characteristic QED value Es = m^c^/eh = 1.32 • 
lO^^V/cm, and thus nonlinear QED vacuum polarization 
effects will become measurable. Some of them have been 
already studied in literature. Among those, harmonic 
generation by an intense laser beam propagating in an 
external magnetic field |j|, and by two colliding beams 
in the vacuum [3. Recently the experimental feasibil- 
ity of the light-by-light stimulated scattering via three 
laser beams has been studied in Ref. [Q|. Pair creation 
by a focused laser pulse in the vacuum was considered in 
Ref. I3, and by two colliding pulses in Ref. [3. Other 
references could be found in the recent review |9|. In 
this letter, we consider the effect of odd harmonics gen- 
eration by a strong focused laser beam in the vacuum. 
Clearly, the probability of this process is smaller com- 
pared to an analogous light-by-light scattering effect in 
a combination of laser beams ^,6] . However it is of great 
importance, since the superhigh laser intensities can be 
achieved only in focused and very short laser pulses, and 
hence the harmonics generation effect will necessarily be 
present at any facility producing pulses with the peak 
field strength of the order of Eg. 
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FIG. 1: a) The exact diagram for generation oil = (2fe + l)-th 
harmonic (a thick line corresponds to the exact electron prop- 
agator); b), c) The diagrams for generation of I = {2k + l)-th 
harmonic in the two leading orders of perturbation theory. 
Dashed lines represent laser photons of frequency lo. 



Feynman diagram depicted in Fig. ^. In the two lead- 
ing orders of perturbation theory the effect is represented 
by diagrams in Fig.^,c. The diagrams without outgoing 
laser photons do not contribute since they would corre- 
spond to creation of an extraneous photon by a group 
of laser photons propagating in one direction. As is 
well known, such photons do not interact. For the fields 
weakly varying at distances of the order of the Compton 
length Z(7, and time intervals Ic/c, nonlinear vacuum po- 
larization effects can be described using the Heisenberg- 
Euler radiative correction to the electromagnetic field La- 
grangian density [lO|, llll Il2l | , which can be represented 
by the following asymptotic expansion 



£' = 



4^2 + 7^ T {8T^ + 13g2 



3607r2^| 



6307r2^| 



, (1) 



a — e^ /he is the fine structure constant, T = (E^— H2)/2 
and t/ = E • H are the invariants of the electromag- 
netic field^. The relation X ^ Ic, where A is the wave 
length of laser radiation, is valid for any of existing or 
projected lasers (with a possible exception of 7, or hard 
X-ray lasers). Thus the Lagrangian density iQJ can be 
evidently applied for description of focused laser fields. 

The Lagrangian density ^ generates effective vac- 
uum polarization P = 9£'/9E and magnetization M = 
—dC /da and the effective vacuum charge and current 
densities 



-VP, j = 



'at 



VxM 



(2) 



Then the modified vacuum Maxwell equations in the 
Lorentz gauge are reduced to the form 



□^^ = 4^J/.(E,H). 



(3) 



The process under consideration is described by the ^ From now on, the natural units h -- 



1 are used. 



The vacuum current (0) contains the fine structure 
constant a. Therefore we will look for a solution to 
Eq. Q (as well as the fields E, H) in the form of per- 
turbative expansion, A^ = A^°'' + A'^ (E = e(°' + E', 
H = H^°-' + H'), where the fields with zero index obey 
the classical vacuum Maxwell equations. For the field in 
the zeroth-order approximation we will use the model of 
a focused laser beam developed in Ref. [13 . This model 
is based on an exact solution of Maxwell equations and 
was successfully employed to explain experimental results 
on scattering of relativistic electrons by a focused laser 
beam 14], and to estimate the number of pairs created 
by a focused laser beam in the vacuum 7, 8]. A similar 
model was used in Ref. [l5J. 

Following Ref. [IJl we describe a focused laser beam 
as a superposition of monochromatic plane waves of the 
same frequency to and with wave vectors lying inside a 
cone with aperture angle 2A. The vector potential of the 
beam propagating along the z axis can be written then 
as 



A("^(rj^,z,i)= f (fkj_a{k_i_)e 



i(k.±r±-\-kzZ — Ljt) 



(4) 



k± <uj 



where 



k = k| 



-kl, k-a(ki)=0. 



The configuration of the beam is determined by the func- 
tion a(kj^) and may vary strongly. However, indepen- 
dently of the form of a(kx), any polarized beam may be 
expressed as a superposition of e- and /i-polarized waves, 
i.e. the waves respectively with the vectors E*^ -' = iujA^ ' 
and H*^"' = V x A^^-* perpendicular to the propagation 
direction of the beam. We will choose the following form 
for functions a.'^'^{'k±) for e- and h- circularly polarized 
beams 



a''^">{ki_) = 



ii;rfci(i±(k)-e)U(ii)(k) 

4\/27rw5A4 



-k±/4uj^A' 



, (5) 



where l_L(k) = (k x z)/k±, l||(k) = (k x 1±)/(jJ and 
e = (x ± iy)/^/2. The resulting expressions for the 
fields can be found in Refs. 0, 113 • They describe a 
focused beam with the focal spot of the characteristic ra- 
dius R — 1/ujA and the diffractive length L — l/wA^. It 
will be important for us later on in this Letter that the 
fields acquire the following structure 



e(o) h(°) 



E^,°^e-^^Fii,x), 



(6) 



where ip — Lu(t ~ z), ^ = r±/R, x — z/L. The quantity 
£"0 is the amplitude of electric (magnetic) field at the 
focal axis of the beam. It is worth noting that the type 
of polarization of the beam is determined by the field 
configuration near the focal axis (see Refs. [ij, [iJI for 
details) . 



For A <C 1 the integral in Q) can be evaluated analyt- 
ically. The formulas for the field invariants J-, Q can be 
found in Ref. [3|. Here we give explicit expressions only 
for the invariants J-, Q in the focal plane z = (0 is the 
polar angle in that plane), 



+6^ cos[2{ujt T (/))] - 28^2 + 4} , 
5"= = 2£;^A2^2g-252 sin[2(wi T </>)], 



(7) 



the upper (lower) sign corresponds to the clockwise- 
(counterclockwise-) polarized waves. The fields in e- 
and /i-polarized waves are related as E^"-* = H*-°^'^, 
ll(o)h = _E(o)e fj^^ -pj^^g^ invariants for the h-wave are 
given by T'^ = —T'^, Q^ = —Q'^- It is clear that in both 
cases they decrease exponentially away from the focal 
spot. Similarly, the invariants decrease fast (in general 
as x~^) away from the focal plane. Therefore, the in- 
duced vacuum polarization charges and currents l(2Jl are 
almost localized inside the focal region with the volume 
Vf ^ R^L ~ l/cj'^A^. For small A, the characteristic 
dimensions of this region are essentially greater than the 
wavelength of the laser field. Hence, space outside of 
the focal region constitutes a wave zone for the radiation 
emitted by the vacuum polarization currents. 

It is convenient to expand the fields E', H' gener- 
ated by vacuum currents in a Fourier series, e. g., H = 
J2i Hj(r)e-*"". It is weh known, see, e.g., Ref. fli|, that 
Fourier components of magnetic strength of the radiation 
field in the wave zone are given by 

T 

■1 iuiir r r , , 

HKr) = ^^^y ^^7 ^'^'^*'"*' -"'■Uxj(r',0, (8) 


n — r/r. Using Eqs. Q we reduce Eq. JSJ after integra- 
tion by parts to the form 



h; 



H^ 



iuilr ( i{^) 



.w 



- f dt' /d3r'e'''^(*'-"'-')/if)(i',r'), 



T 



(9) 



?'(^',r')=(^) l.(n).(E(°)+„xH(°); 

'f)^l.(n).(H(0)-„xE(0)). ^''^ 

The angular distribution of the radiation power at fre- 
quency ujI is given by dli = |Hjpr^ (irik/27r, see Ref. flGj. 
To obtain the rate Pi of photon emission we should di- 
vide this quantity by the photon frequency ojI and then 
integrate it over the spatial angle, 



P; = 



~2^ 



df7„(i^wp+i/.|ii'n 



(11) 



We arrive at the same result computing the Feynman dia- 
gram Fig.^ where the bold line corresponds to the exact 
electron propagator in a locally constant electromagnetic 
field 0. 

Expanding the vacuum current in a Fourier series, we 
represent the probability of photon emission as a sum of 
partial probabilities, each with its own conservation law 
corresponding to conversion of some number n of laser 
photons into n' laser photons (n' < n) and an extraneous 
one with the frequency to' = luj. It immediately follows 
from the Furry theorem that the frequency multiplier I 
can be an odd number only, I = 2k +1. This can be seen 
also from Eq. I|ll|) . since the derivatives of the effective 
Lagrangian (^ with respect to the invariants are even 
functions of electromagnetic field strengths. 

Let us estimate now the rate of photon generation Hll() 
under assumption A <C 1. First we note, that due to 
Eq. ©, hi {t',r') is a function of variables ip' , ^' , (j)' , x' ■ 
Passing on to these variables, we write down Eq. (|10|) in 
the form 



2ir 

hf'> = ^Jd^'Jededcj^'dx'W\v',^'^x') 



f,, , . sin0_^, . ,, .1 -00861 , 

X exp < iLif — zt 4 smcp + it — ^ — x 



(12) 



where 9 is the azimuthal angle of the emitted photon, 
con 6 = (n • z). As it can easily be seen from Eq. |(SJ|, 
the effective values of the variable of integration fc^ in 
the integral Q) are of the order fc^ ^ cjA when A ^ 1. 
Then, it can be shown using Eqs. jsj and the definitions 
of the fields e'°\ H'"^ (gj), that the combinations 



(n) • (: 
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(n) • (: 



e(0) 
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which appear in Eq. (|10|l . are of the order E^ ' H^^ for 
both e- and /i-polarized waves. Both invariants T and Q 
are ~ (A£^q )^ in the focal region, see Eq. ([7||. Then, it 
follows from Eq. Q that the derivatives dC/dJ-^ dC/dQ 
are of the form ag{'d;ip' ,$,' jx') near the focus, where g 
is a dimensionless function, and the dimensionless pa- 
rameter d is defined as ?? = ^-^o l^s- Thus the co- 
efficients ft-j in Eq. © acquire at A ^ 1 the form 



.W - 



h\'' = aVfE^^'i:^^gu{v,^), where v = ^/A. We have 
taken into account the fact that photons are emitted ba- 
sically forward with 6* < A when A ^ 1, as it can be seen 
from Eq. H12|l . Substituting this estimate into Eq. 1)11(1 
we obtain for the number of photons generated per one 
period, Ni = PiT: 



Ni^LjH^a^{Vf)^iEJ,"^AYA^ 



Y, j dvv\gu{vM ^o.(££f ^^im- 



(13) 



The quantity (m/wA)'' has the order of the ratio VfT/l^ 
and is extremely large for optical lasers. For example for 
a laser field with the wavelength A ~ l/im and A = 0.1 
it is ~ 10^^. Therefore the number of generated photons 
will become perceptible already at ^ ^ 1, i.e. at the 
laser field strength E^ essentially less than Es- 

In the framework of the perturbation theory the lead- 
ing contribution to generation of the /-th harmonic [l = 
2fc + 1) arises from the diagram with I + 1 incoming 
and a single outgoing laser photon lines, see Fig. ^3- 
Hence, fi{^) = Ci^^^+^. The coefficients d can be 
computed analytically. For the first few harmonics one 
obtains Cg = 6.380 • 10-^ C5 = 1.048 • IQ-^, and 
Cj = 3.060 • lO^'^. Thus for the number of generated 
photons at z9 = 10^^ we have N3 « 10. For the laser 
field with the wavelength A ~ l/itm and A = 0.1 the 
value ■& = 10~^ corresponds to the peak field strength in 
the focus Eq w O.lEs, or to laser intensity / ~ 10~^/s, 
where Is = (c/47r)£;| = 4.65 • 1029W/cm^ 



It is clear, that at E. 



(0) 



Es and higher, the num- 



ber of generated photons becomes enormous. There- 
fore it makes sense to compare the energy Wt fiowing 
through the focal spot during the wave period T with 
the energy of photons Wrp' generated during the same 
interval of time. The quantity Wt can be estimated as 



Wt ^ {E^' /4:7T)7rR^T, while W^'^' = 'E'^l^i, with Ni 

from Eq. ^. Finally, we have W^'^'' /Wt ^a^J2 ^fiW- 

I 
We see that this ratio remains small at least for i9 close 

to unity. This means that exhaustion of the laser pulse 

due to harmonics generation effect does not occur, as it 

happens in a pair creating laser pulse, see Ref. Q- 

In the region "d ^ 1, where the perturbation expansion 
(^ becomes inapplicable, we have employed the repre- 
sentation of the Lagrangian density by a convergent se- 
ries suggested in Ref. [13 and evaluated the integrals in 
Eqs. (|12|l . 1)11(1 numerically. The results of numerical cal- 
culations are shown in Fig. 12 One can see that, at d close 
to unity, the perturbation theory overstates the number 
of emitted photons. One can see also the tendency to- 
wards equalization of the number of emitted photons in 
neighboring harmonics when -d grows. 

It is worth noting that the probability of the third har- 
monic generation in the considered case is several orders 
larger than that of the second harmonic generation in 
the crossed beam setup p. This can be explained by 
the effect of stimulated emission of the photon. To make 
certain of that, we first notice that both processes are 
described by topologically equivalent Feynman diagrams 
depicted in Fig. lb. The only difference consists in the 
fact that one of the outgoing photons in our case is a 
laser photon, while in the case considered in Ref. pj both 
photons are extraneous. Consider the effect of the third 
harmonic generation by two pairs of laser photons col- 
liding at an angle 2ip, ip < A. At the reference frame 
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FIG. 2: Number of photons generated per wave period by a a) 
e-polarized beam (A^f ) and b) /i-polarized beam (A^;'' ) versus 
parameter I /Is- The sohd curves are the resuhs of numerical 
computation, dashed lines display the results of perturbation 
theory, A = 0.1. 



K' ^ where the collision is head-on, the probability of our 
process per unit time and unit volume can be estimated 
using Eq. (4.38) of Ref. 5]. Since the probability is in- 
variant, we have in the laboratory frame 



dN. 



E-H) 



dVdt 



Es 



dn,i 



(14) 



Here to' and El^ are the frequency and the peak field 
strength of the laser at the X'-frame. The velocity of 
the A''-frame is w ~ cos-0 ~ 1 — V'^/2, and the corre- 
sponding Lorentz factor 7 ~ I/?/;. Hence, oj' — ojip and 
E'q = Efjip. The factor dn^ arises due to the effect of 
stimulated photon emission, and has the meaning of the 
number of laser photons propagating at the interval of 
angles from ^ to ip^ dip relative to the direction of prop- 
agation of the laser beam. It is easy to see using Q that 
dn^ ~ |a(kj_)pdfc^ with fc_L = ujtp, and a(kj^) given by 
Eq. ((SJ. After performing integration over ip in H14|l . and 
multiplying the result by the 4- volume of the focal region 
VfT ~ 1/(0; A)^ we arrive, up to a numerical factor, to 
Eq. 1131) with I = 3. 

In the present letter we have studied the effect of har- 
monics generation by an intense focused laser pulse in 



vacuum. We have shown that for the laser pulse with 
A = 1/xm and A = 0.1 the number of generated photons 
reaches the value of one photon per period at intensity 
J « 5 • lO^^W/cm « 10~^/s. The corresponding peak 
value of the electric field is one order of magnitude less 
than the characteristic QED field Es- This is explained 
by a very large value of the effective focal region as com- 
pared to the characteristic Compton 4- volume, Iq/c. It is 
very important that the rates of photon generation for e- 
and /i-polarized waves are close to each other and reach 
an observable level at z? « 10~^. It follows from Ref. [7| 
that the effect of pair creation by a focused e-polarized 
laser pulse arises at the same values of ■(?, while in an 
/i-polarized wave production of pairs begins noticeably 
later. This means that in an /i-polarized wave we will 
not have the background of secondary photons radiated 
by created particles, and hence, just such waves should 
be employed for observation of the effect of harmonics 
generation by a focused laser pulse in vacuum. 
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